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Enclosure  1 


Quantum  analysis  of  a  micro  cavity- tuned  Bloch  oscillator  for 
tunable  spontaneous  emission  and  absorption  of  terahertz 

radiation 


G.  J.  Iafrate 

Department  of  Electrical  and  Computer  Engineering, 

North  Carolina  State  University,  Raleigh,  North  Carolina  27695-8617 

Abstract 

In  this  study,  we  investigate  the  spontaneous  emission  of  radiation  for  a  Bloch  electron  travers¬ 
ing  a  single  energy  miniband  of  a  superlattice  in  an  external  homogeneous  electric  field  subjected 
simultaneously  the  influence  of  resonant  microcavity  and  dephasing  effects  from  an  internal  in¬ 
homogeneous  electric  field.  It  is  shown  that  the  spontaneous  emission  for  the  cavity-enhanced 
Bloch  electron  probability  amplitude  becomes  damped  and  frequency  shifted  due  to  the  perturb¬ 
ing  inhonrogeneity  when  treated  in  a  long-time ,  time-dependent  perturbation  theory  relative  to  the 
Bloch-accelerated  dynamics  in  the  electrodynamic  radiation  field.  The  frequency  shift  is  shown  to 
be  proportional  to  the  diagonal  matrix  elements  of  the  Hamiltonian  for  the  perturbing  inhorno- 
geneity  with  respect  to  the  instantaneous  Bloch  eigenstates,  and  the  damping  term  is  shown  to 
be  proportional  to  the  off-diagonal  transition  matrix  elements  of  the  perturbing  Hamiltonian  with 
the  instantaneous  eigenstates  summed  to  the  appropriate  final  states  as  determined  in  a  golden- 
rule  like  fashion.  The  resulting  general  theory  is  reduced  for  the  specific  cases  of  an  abrupt  and 
smoothly  varying  potentials  but  emphasis  is  given  to  the  special  case  of  a  comb  of  Slater-Koster 
interface  impurities  with  randomly  distributed  interface  roughness  at  all  lattice  sites.  From  the 
Slater-Koster  case,  the  relaxation  time  approximation  is  developed  where  the  damping  term  is 
considered  to  be  a  constant  and  the  frequency  shift  is  ignored.  Analysis  of  total  power  shows 
that  dephasing  degradation  effects  are  more  than  compensated  for  by  enhancements  derived  by 
microcavity  confinement. 
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I.  INTRODUCTION 


A  theory  of  microcavity  enhanced  spontaneous  emission  (SE)  for  a  Bloch  electron  travers¬ 
ing  a  single  energy  band  accelerating  in  an  external  constant  electric  field  has  been  recently 
examined  by  the  author  and  colleagues.1  The  theoretical  analysis  was  fully  quantum  me¬ 
chanical  in  that  the  quantized  radiation  held  was  described  in  terms  of  the  dominant  TEio 
rectangular  microcavity  waveguide  mode  in  the  Coulomb  gauge;  also  the  instantaneous 
eigenstates  of  the  Bloch  Hamiltonian  were  utilized  as  the  basis  states  in  describing  the 
Bloch  electron  dynamics  to  all  orders  in  the  constant  electric  held.  Analysis  of  the  proba¬ 
bility  amplitudes,  over  integral  multiples  of  the  Bloch  period,  resulted  in  selection  rules  for 
photon  emission  in  both  photon  frequency  and  wave  vector  showing  preferable  transitions  to 
the  Wannier-Stark  ladder  levels;  it  was  shown  that  the  SE  emission  rate  could  be  enhanced 
by  tuning  the  emission  frequency  to  align  with  the  cavity  mode  spectral  density  peak,  thus 
resulting  in  an  output  power  of  several  microwatts  in  the  terahertz  spectral  range  for  a 
GaAs-based  superlattice  imbedded  in  a  microcavity. 

In  this  study,  the  analysis  of  microcavity  enhanced  spontaneous  emission  (MESE)  for  a 
Bloch  electron  accelerating  through  a  single  miniband  of  a  superlattice  structure  was  ex¬ 
tended  to  include  the  additional  interaction  of  the  electron  with  a  perturbing  inhomogeneous 
electric  field  arising  from  impurities  or  interface  roughness  inherent  from  the  superlattice 
material  process.2  The  intent  in  studying  the  effects  of  such  perturbing  inhomogeneities  is 
to  determine  their  role  in  limiting  the  MESE  process  due  to  scattering  influences  that  de¬ 
phase  the  coherency  of  the  Bloch  oscillation,  and  to  establish  a  quantitative  determination 
as  to  how  the  MESE  selection  rules  are  influenced  by  such  perturbations.  For  that  purpose, 
the  Hamiltonian  for  the  Bloch  electron  in  the  quantum  electrodynamic  field  of  interest  is 
developed  (Sec.  II).  The  classical  external  electric  held  is  described  in  the  vector  poten¬ 
tial  gauge,  and  the  quantized  electromagnetic  radiation  held  is  described  by  the  dominant 
microcavity  TE10  rectangular  waveguide  mode  in  the  Coulomb  gauge;  the  general  inhomo¬ 
geneous  electric  held  is  treated  in  the  scalar  potential  gauge.  In  neglecting  the  higher-order 
quantum  field-held  interaction  term,  it  is  shown  that  the  total  Hamiltonian  for  this  problem 
reduces  to  the  sum  of  three  contributions,  the  Hamiltonian  for  the  Bloch  electron  in  the 
classical  external  electric  held  interacting  with  the  inhomogeneity,  the  Hamiltonian  of  the 
free  quantized  electromagnetic  held,  and  the  Hamiltonian  for  the  hrst-order  interaction  be- 
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tween  the  cavity  quantum  field  and  the  accelerated  Bloch  electron.  Then,  the  instantaneous 
eigenstates  of  the  Bloch  Hamiltonian  and  the  states  of  the  free  radiation  field  are  utilized  as 
basis  states  in  describing  the  time  development,  and  in  calculating  the  SE  transition  rates 
of  the  accelerated  Bloch  electrons  under  the  simultaneous  perturbing  action  of  the  quan¬ 
tum  cavity  radiation  field  and  the  inhomogeneous  potential  energy  (Sec.  111).  In  treating 
the  perturbing  inhomogeneity  in  long-time,  time-dependent  perturbation  theory  relative  to 
the  Bloch-accelerated  system  in  the  electrodynamic  radiation  field,  it  is  found  that  the  SE 
amplitude  for  the  cavity-enhanced  Bloch  electron  radiation  becomes  damped  and  frequency 
shifted  in  the  off  diagonal  and  diagonal  matrix  elements  of  the  inhomogeneous  potential 
energy  with  respect  to  the  instantaneous  Bloch  eigenstates.  In  Sec.  IV,  the  general  theoret¬ 
ical  analysis  is  applied  to  the  specific  cases  where  the  inhomogeneous  potential  is  a  slowly 
varying  or,  in  contrast,  an  abruptly  varying  function  of  the  coordinates  as  well  as  a  comb 
of  Slater-Koster  impurities  of  varying  strength  positioned  at  all  the  lattice  cites.  Here  it  is 
found  that  the  frequency  shift  is  a  constant,  and  the  damping  constant  is  slowly  varying 
with  the  Brillouin  zone  vector.  Therefore,  for  purposes  of  showing  trends  with  regard  to 
simple  dephasing  effects,  the  dephasing  analysis  is  carried  out  in  Sec.  V  for  the  case  where 
the  damping  term  is  a  constant  and  the  frequency  shift  is  ignored;  also  a  numerical  estimate 
of  SE  power  is  provided.  In  Sec.  VI,  a  summary  of  the  most  important  results  is  given; 
most  importantly,  our  analysis  of  the  total  power  shows  that  dephasing  degradation  effects 
are  more  than  compensated  for  by  the  enhancements  derived  by  microcavity-based  confine¬ 
ment  engineering.  Section  VII  contains  two  Appendixes,  Appendix  A  and  Appendix  B;  in 
Appendix  A,  we  provide  the  details  of  the  time-dependent  perturbation  theory  analysis  used 
to  calculate  the  probability  amplitude,  and  Appendix  B  provides  the  details  of  calculation 
of  matrix  elements  of  perturbation  Hamiltonian. 

II.  QUANTUM  APPROACH 

The  dynamical  properties  are  considered  for  the  situation  in  which  the  electron  is  confined 
to  a  single  miniband  “no”  of  a  superlattice  with  energy  eno(K)  while  the  effects  of  interband 
coupling  are  ignored.3  Therefore,  the  quantum  dynamics  is  described  by  the  time-dependent 
Schrodinger  equation 

=  (i) 
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where  the  exact  Hamiltoian  H  =  [ p  —  (e/c)A]2/2mo  +  Vc(r)  +  Hr  +  V (r,  t)  can  be  reduced  to 
a  sum  of  the  following  separate  Hamiltonians  H  =  Hq  +  V(r,  t)  +  Hr  +  Hi.  Here,  the  first  two 
terms  represent  the  Hamiltonian,  H0(t)  =  [p  +  pc(t)]2/2mo  +  Vc{r) ,  for  single  electron  in  a 
periodic  crystal  potential,  Vc(r),  interacting  with  a  homogeneous  electric  held  E,  and  V(r ,t) 
is  the  potential  energy  for  the  inhomogeneous  electric  held  due  to  the  impurities  or  interface 
roughness;  also  Hr  is  the  Hamiltonian  for  the  cavity  mode  electromagnetic  radiation  held  (mo 
is  the  free-electron  mass,  c  is  the  velocity  of  light  in  vacuum).  The  total  vector  potential  in 
the  exact  Hamiltonian  consists  of  A  =  Ac  + Ar,  where  Ac  =  — (c/e)pc  describes  the  external 
electric  held  with  pc(f)  =  e  E (t')dt'  =  eEt  for  the  time-independent  homogeneous  electric 
held  turned  on  at  initial  time  t0  =  0;  Ar  describes  the  TE10  cavity  mode  of  the  quantized 
radiation  held  Er  with  the  frequency  uq  given  by 


Ar.y 


I  Anhc2 
LUqeV 


■ in(qxx )  (. 


aqeigzZ  +  a/qe 


t  e-iqzz\ 


where  a'q  and  aq  are  the  photon  boson  creation  and  annihilation  operators,  £  is  the  dielectric 
constant  of  the  medium  filling  the  waveguide  of  the  length  Lz  and  cross  section  Lx  x  Ly, 
V  =  LxLyLz,  and  qx  =  n / Lx.  For  the  chosen  system  geometry,  the  corresponding  electric 
held  Er  is  polarized  in  the  direction  of  the  dc  held,  which  is  assumed  to  be  along  the 
y  axis  (also  the  superlattice  growth  direction).  The  normalization  constant  in  Eq.  (2)  is 
chosen  in  such  a  way  that  the  Hamiltonian  for  the  quantized  radiation  held  has  the  form 
Hr  =  Eg  htuqa^qaq,  where  ujq  =  ujc[1  +  {qz/ Qx)2]1^2  is  the  mode  dispersion  relation,  and 
<^c  =  qxcjs/e  is  the  angular  cutoff  frequency.  The  guided  mode  wavelength  is  written 
as  A  =  \c/[(ujq/ ujc)2  —  l]1/2,  where  Ac  =  2 Lx  is  the  cutoff  wavelength.4  The  Hamiltonian 
Hj(t )  =  —  (e/moc)Ar-[p+pc(t)],  for  the  hrst-order  interaction  between  the  quantum  held  and 
the  Bloch  electron,  couples  both  subsystems  H0  and  Hr,  and  causes  transitions  between  the 
accelerated  Bloch  electron  states  through  photon  absorption  and  emission.  Then,  starting 
with  the  reduced  Hamiltonian  H  =  H0  +  F(r,  t)  +  Hr  +  Hj,  use  is  made  of  hrst-order,  early- 
time  time-dependent  perturbation  theory5  to  calculate  SE  transitions  probabilities  between 
states  of  H0  +  Hr  while  regarding  Hj{t )  ~  Ar-  [p  +  pc(t)]  as  a  perturbation,  and  at  the  same 
time,  use  is  made  of  a  long-time  perturbation  theory  analysis6  to  calculate  the  relaxation 
influence  of  V (r,  t )  on  the  SE  transition  probabilities.  This  time  scaling  is  justified  since  SE 
rates  are  orders  of  magnitude  faster  than  relaxation  times  due  to  V{r ,t). 
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III.  QUANTUM  DYNAMICS  BASED  ON  INSTANTANEOUS  EIGENSTATES 

The  solution  for  |Tno(f)}  in  Eq.  (1)  can  be  represented  in  terms  of  the  eigenstates  of  basis 
states  |^nok (t),  {nqj})  =  |Vu0k(t))  |{nqj})  of  the  unperturbed  Hamiltonian  H0  +  Hr  as 

i«™«>  =  E  E  A*.  j}(k>  (t)i 

k  {nq5j } 

^  [£n0(k(t'))  +  j  >  (3) 

where  the  summation  over  k  is  carried  out  over  the  entire  Brillouin  zone,  and  {nqj}  is 
specified  over  all  possible  combinations  of  photon  occupation  number  nqj-  with  photon  wave 
vectors  q  and  polarization  j  =  1,  2.  The  instantaneous  eigenstates  of  H0  are  given3  by 
'0nok(t)(r,  t)  =  Q_1/2e*K'rMnok(t)(r,  £) ,  where  wnok(t)(r>  £)  is  the  periodic  part  of  the  Bloch 
function,  k(f)  =  K  +  p c(t)/h,  and  the  values  of  the  electron  wave  vector  K  are  determined 
by  the  periodic  boundary  conditions  of  the  periodic  crystal  of  volume  Q. 

For  the  case  of  one-photon  SE,  which  assumes  that  initially  no  photons  are  present  in  the 
field,  the  probability  amplitude  in  the  wave  function  of  Eq.  (3)  satisfies  the  initial  condition 
A{nq  }(k,  t0)  =  {A)q  ,,0}  (5k, K0  at  time  t  —  t0  when  the  electric  field  is  turned  on.  Here,  K0 
and  nq  •  =  0  are  the  initial  values  of  K  and  nqj- .  Thus,  the  instantaneous  eigenstates  of  H0 
satisfy  the  equation3 

[p  +  Pc(t)]2  +  K(r)  j  V’nok(t)  =  £no(k(t))^n0k (t)  ,  (4) 

and  the  basis  states  for  Hr  are  given  by  the  well-known  free  photon  field  equation 

Hr\{nc\.j})  —  'y  ^  flUJqin  qj- 1  {n.qj  } )  ,  (5) 

q  ,j 

where  {nqj})  is  a  simple  product  of  all  possible  combinations  of  photon  number  states, 
nqj,  with  a  given  wave  vector  q  and  polarization  £qj . 

The  appropriate  time-dependent  equations  of  motion  for  the  A{„qi}(k,  t)  coefficients  ex¬ 
pressed  in  Eq.  (3)  can  be  found  in  Appendix  A;  these  equations  relate  the  time  rate  of 
change  of  A{n  }(k,  t)  to  the  basis-dependent  matrix  elements  of  Hj  and  V(r,t)  through  a 
self-consistent  set  of  equations.  In  applying  early-time,  first-order  perturbation  theory  to 
the  ”  Hi”  coefficients  of  the  set  of  equations,  and  applying  a  long-time,  Wigner-Weiskopf-likc 
approximation  to  the  ”V”  coefficients  of  the  set  of  equations,  we  obtain  (see  Appendix  A)  a 
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closed  form  inhomogeneous  equation  for  the  one-photon  SE  amplitude  [here  |nqj}  — >  n  —  1] 
at  k0(t)  =  K0  +  p c(t)/h  as 


1  l 

An(k0,t)  — /k0ko(t)^4n(ko,  t)  + 


/  dt>  /k()k'  W/k0k'  (^)  Ai(k0,  t') 
'  °  kVko 


(6) 


k'  /k0 


Here, 

/kok'W  =  ^k0k'  exp  /  KoMO]  -  £n0[k'(t,)]}^/}  ,  (7) 

where  Vkok/  =  (VVioMth  ^V’nok'(t))  are  the  appropriate  matrix  elements  of  the  perturbing 
inhomogeneity  U(r,f),  and 


Hj(k0,t)  =  D(qx/q)1/2  /  df'uJkoO')  -  qs] 


x  exp  —  —  l  {£n0  [ko(H)]  £n0  [k0(H)  -  q.s]  -  fuVqjdti  )  , 


(8) 


the  probability  amplitude  for  the  microcavity-based  SE  alone.1  Also,  D  =  —i^/nca/uJeeV, 
a  =  e2/ he  is  the  fine  structure  constant,  qs  =  {±g.T,  0,  qz}  with  ”  +  ”  for  s  —  1  and  ”  —  ”  for 
s  =  2,  vy(k(t ))  =  (1  /h')Vj£y£no (K ) |k(i) j  the  y  component  of  Bloch  velocity  in  the  band,  and 

q  =  ( ql  +  ql)1' 2- 


IV.  ANALYSIS  OF  DEPHASING  FOR  PARTICULAR  MODELS  OF  INHOMO¬ 
GENEITY 

The  solution  for  An(k0,f)  in  Eq.  (6)  depends  explicitly  on  the  potential  energy  V{r ,t)  , 
through  /k||k'(t)  as  noted  in  Eq.  (7).  In  general,  the  time  dependence  of  /k(lk'(t),  expressed 
through  the  time  dependent  matrix  elements  of  V (r,  t )  and  through  the  time  dependence 
of  accelerated  instantaneous  energy  eigenstates  in  the  phase  of  /k()k'(t),  makes  finding  the 
solutions  for  An(k0lt)  in  Eq.  (6)  formidable.  However,  for  the  type  of  potential  energy 
function  of  interest  in  this  problem,  that  is,  one  where  V(r ,t)  =  U(r),  a  function  of  position 
alone,  and  where  V (r)  is  a  comb  of  abruptly  changing,  Slater-Koster  like,  interface  impurities 
with  randomly  distributed  interface  roughness,  it  can  be  shown'  that  /kk'  (t)  is  expressed  as 

/kk'(t)  =  IW  e-sWK)-AK')]  ei[£-„0(K±)-^0(K 'x)k  (9) 
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It  is  important  to  note  that  EKK',  77  ( K),  and  £no(K_i_)  are  all  time  independent.  EKK'  are 
the  matrix  elements  of  the  interface  inhomogeneities;  77 ( K)  and  £no(K^)  are  band  structure 
dependent  only ,  and  are  given  by 

v[K)  =  2  6no (aly > Kj- )  gin (Kyaly) ,  (10) 

j„=i 

with 

£„„K,  Kj.)  =  4-  53 Kilt-*'"1' ,  (11) 

y  K, 

and 

£no(K_|_)  =  ——  y  '  £no  (Ky,  Kj_)  ,  (12) 

y  Ky 

the  average  value  of  £no(K)  along  the  Ky  direction  in  the  Brillouin  zone,  independent  on 
Ky\  Ny  is  the  number  of  lattice  sites  along  the  y  axis.  For  nearest-neighbor  tight-binding 
considerations,  only  the  term  with  ly  =  1  need  be  retained  in  Eqs.  (10)  and  (11).  Then, 
in  utilizing  /kk/(f)  from  Eq.  (9)  in  Eq.  (6),  and  anticipating  the  ensemble  average  over  the 
randomly  distributed  interface  roughness  coordinates,  which  renders  the  ensemble  average 
over  the  /kk/(f)  term  on  the  right-hand  side  of  Eq.  (6)  equal  to  zero,  then  the  ensemble 
averaged  equation  for  An(k0,f)  reduces  to 


X(ko,f)-lykoiwA,(k„,t)  +  i  v  |vtotf 

kVko 

ft  _  , 

x  /  es^"o(-Ko-L)_^no(K,-L)](t_t  )j4n(ko,f/)  =  A°(ko,£) .  (13) 

Jo 

This  is  an  inhomogeneous  differential  equation  in  time  for  An(k(),  t),  where  the  inhomogeneity 
on  the  right-hand  side  is  the  time  derivative  of  the  SE  for  the  cavity-based  result  alone  [see 
Eq.  (8)].  This  equation  can  be  readily  solved  by  the  Laplace  transform  method  in  the 
long-time  limit  to  obtain 

An(ko,  t)  =  f  dt  iJ(ko,  t')  e~iAu<*>)*  ,  (14) 

Jo 


where 


and 


Acu(ko)  -  4oko  +  r 


IH 


k(,k 


2vr 


^  •  ^no(K-OJ-)  (K  _L y 


e|7m[,?(K0)-r,(K  )] 


kVko 


(15) 


r(k»)  =  lT  E  lUokf  e*,’"l”(K“>-’'<K)15[£„(Kox)-£„0(K'±)] 


k'A» 


(16) 


This  result  for  An(ko,  t)  in  Eq.  (14)  shows  that  the  SE  for  the  cavity-enhanced  Bloch  electron 
probability  amplitude  becomes  damped  by  T(k)  in  Eq.  (16)  and  frequency  shifted  by  Acn(k) 
given  in  Eq.  (15)  dne  to  the  perturbing  influence  of  the  assumed  interface  inhomogeneity 
and  the  spatially  localizing  effect  of  the  constant  electric  field  from  r)( K)  defined  in  Eq.  (10). 
It  is  interesting  to  note  that  the  result  obtained  in  Eq.  (14)  is  strikingly  reminiscent  of  the 
result  obtained  for  the  solution  of  the  kinetic  Boltzmann  equation8  for  the  carrier  distribu¬ 
tion  function  in  the  relaxation  time  approximation.  Therefore,  in  an  approach  similar  to 
analyzing  the  complications  of  the  Boltzmann  theory,  we  first  consider  the  model,  noted  as 
the  relaxation  time  approximation,  where  A uj  =  0  and  T(k)  is  a  constant,  independent  of 
k,  to  analyze,  in  the  simplest  heuristic  approximation,  the  degrading  effects  of  dephasing. 


V.  RELAXATION  TIME  APPROXIMATION 


In  this  section,  for  purposes  of  showing  heuristic  trends  with  regard  to  simple  dephasing 
effects,  the  dephasing  analysis  is  carried  out  for  the  case  where  the  damping  term  is  a 
constant  and  the  frequency  shift  is  ignored.  Then,  it  follows  from  Eqs.  (14)-(16)  that  the 
one-photon  SE  probability  amplitude  [i.e.,  {nqj}  — >  n  —  1]  can  be  presented  as 

Aq(k0,t)  =  D(qx/q)1/2  f  dt' vy( k0  -  qs)e“(t'“to)/T 

Jto 


xexp 


£no(ko  — 


qs)  -  hujq\dt 


(17) 


where  r  has  the  meaning  of  the  characteristic  mean  dephasing  time.  Equation  (17)  can  be 
formally  rewritten  as 


Aq(k0,t)  =  D(qx/q)1/2  /  dtf  vy( k0  -  qs) 


' 10 


xexp  <;  --  /  [eno(ko)  -  £n0(ko  -  qs)  -  tujq]dti  , 
11  Jt0 


(18) 


with  a  renormalized  complex  photon  '’energy”  hooq  =  huq-\-ihujT ,  where  cur  =  1/r.  Thus  the 
emission  process  results  in  the  SE  probability 


q  s= 1,2 


(19) 
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A.  Selection  rules 


In  evaluating  Ag(k0,t),  we  take  into  account  that  the  external  dc  field,  E,  is  along  the  y 
axis;  then,  it  follows  that  k0y(t)  =  K0y  +  eE{t  —  t0)/h.  In  taking  advantage  of  the  periodic 
properties  of  the  terms  in  Eq.  (17),  Ag(k0,t)  is  evaluated  in  clocked  integral  multiples  of 
the  Bloch  period,  so  that  t  =  Ntb,  where  Tb  =  2tt /ub,  the  time  to  traverse  one  period  of 
the  Brillouin  zone.  The  integral  in  Eq.  (17)  over  time  can  be  replaced  by  an  integral  over 
koy  through  the  substitution  dt  =  ( h/eE)dkoy .  Then  the  probability  amplitude,  at  integral 
multiples  of  the  Bloch  period,  can  be  expressed  through  that  over  the  single  Bloch  period, 
Tb-3’9  Thus  we  obtain 


Aq(k0,  Ntb) 


1  —  exp(—iN/3q ) 
1  -  exp(-i(3q) 


Aq(k0,rs) , 


(20) 


where  the  complex  parameter  (3q  is  given  by 

q  1  rKoy+Gy 

Pq  =  2vr—  +  —  /  dk0y  [£no(k0)  -  £no(k0  -  qs)]  ,  (21) 

u)B  eE  JKoy 

and  Gy  =  27r /a,  the  y  component  of  the  SL  reciprocal-lattice  vector.  From  Eq.  (20)  the 
corresponding  squared  probability  amplitudes  are  expressed  as 


l4,(ko,ArJs)|2  =  V {w,;N,t)  |A,(k0,  tb)|2  , 


(22) 


where  the  ”  transfer”  function  r/ (ojq ;  N,  r)  is  given  by 


V{Vq,  N,  T ) 


sh2{nN /tub)  +  sin\nNujq/ujB)  c^(N-i)/Tion 
sh?{ii /tub)  +  sin2(nuq/uB) 


(23) 


From  Eqs.  (22)  and  (23),  it  is  seen  that  for  the  limiting  case  of  ttN/tub  -C  1,  i.e.,  when 
the  electron  scattering  is  not  essential  even  for  many  Bloch  oscillations  ( N )  in  the  time  scale 
of  the  scattering  time  r,  the  quantity  \Aq(ko,  Ntb)\2  will  reach  its  maximum  growth  value 
when  j3q  =  2n (m  +  5),  where  f3q  =  Re(/3q),  m  is  an  integer  and  5  — ■>  0;  for  this  limit,  the 
function  r](uq,  N,  r)  is  reduced  to  the  function  r)(uq,  N)  previously  obtained  in  the  absence 
of  scattering:9  rj(uq,N)  =  sm2(N/3q/2)/  sin2(/3g/2)  — >  N2,  i.e.,  it  becomes  sharply  peaked 
at  the  resonances  with  increasing  N.  It  is  clear  that  this  condition  for  maximum  growth 
establishes  the  selection  rule3,9  for  both  the  photon  emission  frequency,  uq,  and  the  key  wave 
vector  component,  qz.  Indeed,  from  the  condition  /3q  =  2n m,  it  follows  from  Eq.  (21)  in  the 
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radiative  long- wavelength  limit  ( qa  <C  1)  that  one  generally10  has 


U) 


g 


mixiBi 


(24) 


this  shows  the  “Stark  ladder ”  resonance  frequency  condition,  along  with  the  sustaining 
wave  vector  cavity  resonance  conditions.  Thus,  the  modes  that  radiate  with  the  highest 
probability  correspond  to  the  fundamental  Bloch  frequency  and  its  harmonics.  These  quan¬ 
tization  conditions  are  obtained  without  requiring  any  assumptions  concerning  the  existence 
of  Wannier-Stark  energy  states.  In  Figs.  1  and  2,  we  show  the  modification  of  the  behavior 
of  r}(ouq,  N,  t )  with  increasing  the  intensity  of  scattering  (decreasing  r)  calculated  for  reason¬ 
able  scattering  ( Ntb/t  =  1)  and  strong  scattering  ( Ntb/t  =  4),  respectively,  with  N  =  10. 
It  is  seen  that  the  peak  positions  are  negligibly  affected,  although  the  peak  spectral  width 
becomes  broadening  and  the  peak  value  is  reduced  with  decreasing  ratio  t/tb,  as  shown  in 
Figs.  3  and  4.  From  this  result,  it  is  noted  that  the  maximum  growth  condition  for  the  prob¬ 
ability  amplitudes  in  Eqs.  (22)  and  (23)  still  establishes  the  selection  rules  for  spontaneous 
emission  given  in  Eq.  (24)  despite  the  increasing  intensity  of  scattering  ( tub  ^  1). 

It  is  worth  to  emphasize  that  in  what  follows  the  equation  (22)  should  be  considered 
simultaneously  with  the  equation  (24),  i.e.,  both  established  resonance  conditions  for  the 
emission  photon  frequency  and  wave  vector  should  be  taken  into  account  explicitly  in  eval¬ 
uating  the  spontaneous  emission  probability  spectra  given  by  |A9(k0,  tb)|2. 


B.  General  expression  for  the  total  SE  probability 


The  total  spontaneous  emission  probability  is  evaluated  at  time  t  =  Ntb ,  P*  =  P*(Ntb), 
by  substituting  \Aq(k0,  Ntb)\2  from  Eq.  (22)  into  Eq.  (19).  The  sum  over  q  in  Eq.  (19)  has 
been  replaced  by  an  integral  over  q,  taking  into  account  the  TEi0  mode  density  of  states 
and  polarization  such  that  ^  (•  •  •)  — >  (Lz/n)  f  dq( ■  ■  •)/[!  —  (qx/d)2]1^2-  Thus  we  can  write 


ps  _  _ ~ 

6  7T 


T  . 2 


z,. 

5=1 


ldqrj(u>q,  N,  r ) 


\Aq(k0,TB)\- 

t1  -  idx/dY 


1 1/2 


(25) 


The  integral  can  be  evaluated  by  using  the  property  of  the  integrand  which  contains  a  sharply 
peaked,  symmetric  function  of  q,  r)(u jq;N,r),  at  q  =  qm  =  (q2  +  q2zm)1^2  [see  Eq.  (24)].  We 
note  that  with  decreasing  the  scattering  time  r,  the  maximum  values  of  this  function  are 
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suppressed  relatively  its  peak  value,  N 2,  corresponding  to  an  infinite  r  (Figs.  3  and  4). 
However,  its  sharp  and  oscillatory  behavior  still  preserved  for  the  scattering  time  not  much 
less  than  the  inverse  Bloch  frequency  uB,  i.e. ,  for  ruB  ^  1  (Figs.  1  and  2).  In  what  follows  we 
assume  such  values  of  r  and  uB  that  satisfy  this  inequality.  Thus,  at  every  node  defined  by 
the  resonance  conditions,  the  slowly  varying  function  of  q  in  the  integrand  can  be  replaced 
by  its  value  evaluated  at  q  =  qm,  and  then  removed  from  the  integral  over  q:  after  that,  we 
obtain 


(26) 


Here  lmax  follows  from  qmax  =  lmax(^B/^c)qx,  and  determines  the  upper  limit  in  the  sum 
over  higher  Bloch  oscillation  harmonics.  The  integral  In  (tub)  is  determined  by 


In  (ruB)  =  e-dv-d/-B_ 

71 


T/<2  ,  sh2(N/2ruB)  +  sin2(Nx) 

nqr - 

sh2(l/2ruB)  +  sin2(x) 


(27) 


The  calculation  of  Pes  in  Eq.  (26)  requires  the  use  of  Aq(k0,TB)  in  Eq.  (17),  evaluated 
at  the  maximum  growth  conditions  of  Eq.  (24),  that  is  when  Uu>q  =  rnhuj B  and  q  =  qm.  In 
addition,  the  dependence  upon  q  in  Eq.  (17)  is  made  explicit  by  invoking  the  assumption 
of  photon  long- wavelength  limit,  which  is  valid  for  all  periodic  potentials  of  interest,  even 
superlattices,  where  q  -C  vr/a.  Thus,  we  find  that 


1/2 


Aqi( k0,  tb)  =  f  —  )  DIi  exp(iK0yauqi / u>B ) , 

uB  \Qi  / 


(28) 


where 


/)  =  -*-/  d'dk[vy('dk)exp('dk/2'KTuB)]exp(-iWk)  (29) 

is  the  Mil  Fourier  component  of  the  function  vy('()k)  exp(i)k/2nTUB),  t^k  =  koya,  and  ujqi  = 
Iub+i/t.  Since  the  electron  velocity  component  vy(flk)  depends  on  details  of  the  superlattice 
miniband  structure,  such  dependence  comes  explicitly  into  the  SE  probability  amplitude 
through  the  integral  Ii  given  by  Eq.  (29);  thus 


|H9;(k0,TB)|2  =  Qx  \DIi\2  exp(-K0ya/nTuB)  •  (30) 

Qi  (cub/ 2tt)z 


C.  Analysis  of  the  total  SE  probability  for  specific  energy  miniband 

The  analysis  for  spontaneous  emission  and  radiation  characteristics  is  now  developed 
by  considering  a  quite  general  form  of  the  electron  energy  miniband  dispersion  relation 
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expressed  as 


OO 


£no(K)  =£no(0)  +  ^A//  si 


laK „ 


sin- 


+  £±(  Kl), 


i'=i 


(31) 


where  £no(0)  is  the  miniband  edge,  A is  the  width  of  the  Z  -th  miniband  harmonic  of 
the  superlattice,  and  £_l(Kj_)  is  the  contribution  from  the  perpendicular  components  of  the 
band.  Such  form  of  the  energy  band  dispersion  in  the  superlattice  growth  direction  generally 
includes  long  range  coupling  over  the  neighboring  QWs  with  a  relative  strength  measured  by 
the  specific  value  of  the  ratio  Aj'+1/Az'  <  1,  which  is  strongly  dependent  upon  the  extent  of 
wave  function  overlap.  In  particular,  the  well-known  case  of  nearest-neighbor  tight-binding 
(NNTB)  approximation  corresponds  with  purely  harmonic  energy  dispersion  where  only  the 
single  term  (/  =1)  is  considered  significant,  so  that  next  nearest  neighbor  and  longer  range 
QW  wave  function  overlaps  are  assumed  to  be  negligibly  small.  The  electron  group  velocity 
in  the  general  miniband  of  Eq.  (31),  for  the  given  Ky  in  the  y  direction,  is  then  given  by 
Vy(Ky)  =  (1  /h)(d£no(Ky) /dKy)  =  Y1T=\  v i'  sin ( l' aKz ) ,  where  ty  =  l'aAl'/2h,  the  maximum 
velocity  associated  with  the  Z  -th  miniband  of  band  width,  A ^ .  Substituting  the  expression 
for  Vy(Ky)  in  Eq.  (29),  one  can  find 


h  — 


1 

2ni 


-  '  >i,  I  Jr;  ..  , 


(32) 


where 


vt>  f  cos[k{1'  —  /)] 


OO  f 

vi  {i{i  - 

i'= i 


cos[7r(//  +  /)]  1 

/)  +  1/27TTZ/B  i(l'  +  /)  —  1/27TTZ/B  J 


Then  the  probability  amplitude  in  Eq.  (30)  can  be  obtained  as 


\Aqi(k0,TB){2  =  ^DStfsh2 
Qi^b 


2  TV 


B 


exp 


I\.Q  yCl  \ 


VT TUB  J 


The  total  spontaneous  emission  probability  is  calculated  from  Eqs.  (26)  and  (34)  as 


P:  =  2aIN(TVB)s'*^Vj-Ash*( 

J-jy  C  /(  y 


2  tvb 


exp 


jy  \  Imax 

J^OyCi  \  _ 

ktvb)  ^  l[  1  - 


IS, 


(33) 


(34) 


.  (35) 


Then,  it  follows  that  for  the  pure  harmonic  miniband  (l' 
reduced  to  the  expression 


1)  the  integral  in  Eq.  (32)  is 


_ sh(l/2ruB) _ 

1  [(1/2 TVs)2  —  vr 2{l2  —  1)]  —  ml /tub  ’ 


(36) 
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thus  to  obtain  from  Eq.  (34)  for  the  probability  amplitude 


/,  x,2  _  2\p\2ucvl  47r2s/r2(l/2rz/B)  exp(—K0ya/ tttub) 

1  0,TB)l  _7F|  1  (ttI/tUb)2  +  [(1/2tUB)2  —  7T2(/2  —  l)]2  ’ 

here  account  has  been  taken  for  the  wave  vector  qi  =  qxluB/cuc,  and  V\  =  aAi/2 h.  In 
particular,  approaching  r  to  infinity,  we  obtain  the  expressions  characteristic  for  the  limiting 
case  when  the  scattering  is  totaly  ignored,  i.e. ,  Ii  =  —i{yi/2)5n  and 


\Aqi^TB)\2  =  AD 


U>B 


Thus,  we  arrive  to  the  conclusion  that  for  the  electron  dynamics  in  a  purely  harmonic 
miniband  when  only  single  Fourier  component  of  the  electron  velocity  vy(i}k)  with  /  =  1  is 
nonzero  and  in  the  absence  of  scattering  (r  — >  oo),  all  probability  amplitudes  |  (k0,  r-^)  | J 
are  zero  except  the  one,  l  =  1,  corresponding  to  the  fundamental  Bloch  frequency  according 
to  the  selection  rule  established  by  Eq.  (24).  However,  with  finite  r  values  the  higher  Bloch 
harmonics  generation  becomes  possible  even  in  the  pure  harmonic  superlattice  miniband. 
It  is  seen  that  with  increasing  r,  the  effect  of  higher  Bloch  harmonics  generation  in  the 
pure  harmonic  miniband  is  suppressed,  and  for  r  — >  oo  it  is  totaly  disappeared.  In  noting 
the  contribution  from  higher  harmonics  of  Bloch  frequency  into  the  total  SE  probability, 
note  that  the  NNTB  approximation  in  Eq.  (31)  is  obtained  by  letting  =  Ai 5U>  so  that 
Vf  =  v\ 8lt'.  The  occurrence  of  the  Kronecker  symbol,  8U> ,  allows  the  contribution  of  l'  =  1 
term  only  whereas  all  the  other  terms  are  equal  to  zero,  thereby  limiting,  within  the  NNTB 
approximation  and  without  of  scattering,  the  generation  to  the  fundamental  Bloch  harmonic. 

From  equations  (26)  and  (37),  the  total  spontaneous  emission  probability  becomes 


_o  1/2  Lx  v\  u2c  4n2sh2(l/2TvB)exp(-K0ya/nTvB) 

e  «  n{tub)  e  ^  ^  ^  2_,  l{1  (£Jc/Z£Jb)2]i/2  {(nl/TUB) 2  +  [(l/2r^)2  -  vr2(Z2  -  l)]2}  ' 

(39) 


D.  Total  SE  power  estimate 


It  follows  from  Eqs.  (39)  that  spontaneous  emission  probability  corresponding  to  the 
fundamental  Bloch  harmonic  (/  =  1)  is  given  by 

r>s(] no  2  nT  /  „  x  v\  u2c  sK2(1/2tvb)  exp{-Koya/ tttub) 

i e  (  w  CX.1^j\TUqJ  6  2  2  M  /  \211/2IY  /  \2  i  /l  /o  \4l  *  \  / 

Ly  c2  ujb  [1  —  u>c/ujB)-j\1'2[(n/Ti,B)2  +  (1/2tvb)  q 
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In  noting  that  the  SE  probability  of  Bloch  radiation  in  the  absence  of  scattering  is  given1 


by  the  expression 


T  ?,2  pl/2,  ,2 

Pse(l  =  1, r  -  oo)  =  2aN  -Z  4  -  , 

Ly  c2  u2B(l-u2c/u2B)V2 


we  can  compare  both  the  probabilities  for  SE  at  fundamental  Bloch  frequency  into  the 
waveguide  TE10  mode  with  and  without  the  scattering  taking  into  account  analyzing  the 


=  Pse{l  =  1,t)  =  1_  sh2{l/2TVB) 

Pg  (l  —  1,t  — >•  oo)  N  N  B  (1/2tz/b)2[1  +  (l /An tub)2}  ' 


here  use  has  been  made  of  Eq.  (40)  for  the  P/(l  —  A ,r)  with  zero  value  of  the  electron  initial 
wave  vector  component  in  the  direction  of  the  constant  electric  field.  From  this  equation, 
we  obtain  for  the  limit  of  r  — ■>  oo,  with  N  finite,  the  expected  result  P/e  — >  1.  Also,  for 
this  limit,  the  integral  (27)  becomes  a  linear  function  of  N ,  In(tub  —■ ►  oo )  =  N ,  and  the 
result  of  Eq.  (40)  goes  to  that  of  Eq.  (41).  For  finite  values  of  the  dephasing  time  r,  it 
is  seen  from  Eqs.  (40)  and  (41)  that  the  linear  dependence  of  P/  on  N  in  Eqs.  (41)  is 
replaced  by  the  factor  In  (tub)  which  is  a  slower  increasing  function  of  N.  The  dependence 
of  In  (tub)  on  N  is  shown  in  Fig.  5  for  several  values  of  the  parameter  tub ■  In  Fig. 
6,  we  show  the  dependence  of  the  relative  SE  probability  PJe  [Eq.(42)j  calculated  as  a 
function  of  number  of  Bloch  oscillations  N  at  different  values  of  dephasing  parameter  tub- 
It  is  seen  that  the  degrading  effect  of  dephasing  on  the  total  SE  probability  within  the 
relaxation  time  approximation,  P is  that  it  becomes  clamped  with  the  extent  of  damping 
dependent  on  the  parameter  tub-  For  numerical  estimation  of  the  damping,  we  assume 
that  tub  =  10  and  N  =  10  which  corresponds  to  the  case  of  reasonable  scattering  with  the 
characteristic  parameter  tub/N  =  1.  Then  using  the  data  of  Fig.  6  (upper  curve),  we  obtain 
P/e  ~  0.5,  that  is  the  total  SE  probability,  and  thus  the  generated  power  of  spontaneous 
emission,  is  damped  approximately  by  a  factor  of  2  due  to  dephasing  in  comparison  with 
that  in  the  absence  of  scattering.  Then,  assuming  the  previously  analyzed  conditions1  for 
the  superlattice  parameters,  the  applied  electric  field,  and  the  waveguide  parameters,  we  can 


use  the  estimated  SE  power  output  in  the  absence  of  scattering,  W  ~  3  /ilT,  which  results  in 
the  estimation  of  the  damped  SE  power  generated  under  dephasing  effects  as  W  ~  1.5  jiW . 
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VI.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 


•  The  quantum  electron  dynamics  and  spontaneous  emission  of  radiation  for  a  Bloch 
electron  traversing  a  single  energy  miniband  of  a  superlattice  accelerating  under  the 
influence  of  superimposed  constant  external  and  inhomogeneous  internal  electric  fields 
while  radiating  into  a  microcavity  have  been  analyzed.  The  analysis  is  based  on  the 
use  of  instantaneous  eigenstates  of  the  Bloch  Hamiltonian  good  to  all  orders  of  the  dc 
field,  and  to  first-order  perturbation  theory  in  the  quantized  radiation  field. 

•  It  is  shown  that  the  spontaneous  emission  for  the  cavity-enhanced  Bloch  electron 
probability  amplitude  becomes  damped  and  frequency  shifted  due  to  the  perturbing 
inhomogeneity  when  treated  in  a  long-time,  time-dependent  perturbation  theory  rela¬ 
tive  to  the  Bloch-accelerated  dynamics  in  the  electrodynamic  radiation  field. 

•  It  is  found  the  frequency  shift  is  proportional  to  the  diagonal  matrix  elements  of  the 
Hamiltonian  for  the  perturbing  inhomogeneity  with  respect  to  the  instantaneous  Bloch 
eigenstates,  and  the  damping  term  is  proportional  to  the  off-diagonal  transition  matrix 
elements  of  the  perturbing  Hamiltonian  with  the  instantaneous  eigenstates  summed 
to  the  appropriate  final  states  as  determined  in  a  golden-rule  like  fashion. 

•  The  analysis  for  the  special  cases  (slowly  or  abruptly  varying  inhomogeneities  as  well 
as  a  comb  of  Slater-Koster  impurities  with  varying  strength  positioned  at  all  lattice 
sites)  showed  that  the  frequency  shift  is  a  constant,  and  the  damping  term  is  slowly 
varying  with  Brillouin  zone  wavevector. 

•  The  general  analysis  allowed  to  justify  the  relaxation  time  approximation  where  the 
damping  term  is  considered  to  be  a  constant  and  the  frequency  shift  is  ignored;  in 
this  case,  it  is  shown  that  the  selection  rules  for  spontaneous  emission  are  negligibly 
affected,  although  the  total  probability  for  emission  shows  interesting  damping  and 
mixing  characteristics. 

•  The  spontaneous  emission  rate  in  the  terahertz  frequency  range  has  been  analyzed  as 
a  function  of  the  Bloch  frequency,  the  external  electric  field,  number  of  oscillations, 
and  dephasing  parameter.  The  analysis  of  the  total  power  shows  that  the  dephasing 
degradation  effects  are  more  than  compensated  for  by  the  enhancements  derived  by 
microcavity-based  confinement  engineering. 


16 


VII.  APPENDIXES 


APPENDIX  A:  TIME-DEPENDENT  DOUBLE  PERTURBATION  THEORY  AP¬ 
PROACH 

The  equation  for  the  probability  amplitudes  A{n  .}(k,  t)  in  Eq.  (3)  is  given  by 

dA["-f{k't) E  A{<  ,(k', ()<{»„■}, v-„0k(,)l(ff,  +  r)ivnst.(t), {«„,}} 


k  K.J 


xexp  <j  /  [£no(k'(t'))  -  £no(k(t'))  +  J  ~  n^Tiu^dt'  }>  .  (Al) 


'to 


qj 


We  assume  that  at  initial  time,  to,  the  system  is  in  one  of  the  eigenstates  of  Hamilto¬ 
nian  H0  +  Hr  with  wave  function  ]V;n0K0,  {nqj})>  corresponding  to  the  Bloch  electron  in 
a  single  band  “no”  with  the  wave  vector  K0,  i.  e.,  V;n0K0  =  (I/O1/2)  e*K°'rwnoKo,  and 
with  the  initial  distribution  of  photon  numbers  in  the  radiation  field  {nW}).  Substituting 
A{„qi}(k,t)  =  |(k ,t)  +  }(k,t)  +  ...  into  Eq.  (Al),  and  taking  into  account  the 

initial  condition  A{„q ^(k,  f0)  =  {Aiq;,:n°  .  }5k(q,k0(t)>  one  obtains  to  the  zeroth  and  first  order 
in  Hi  for  Aj°^_.j(k,  t)  and  A|^^(k,  t),  respectively, 


{^nq,j,nU}^k(t),k0(t)  7 


(A2) 


and 


{nq,j}^  ’  ^  k'e^(-T  [  l£n0(k' (t'))  -  £no{k(t'))]dt' 


dt 


h 


=  ^({n^}7^noHt)\Hl\^noU'0(ty 


xexp  <j  -l-  [£no(k0(0)  “  £n0(k(t'))  +  ^(n-qj  -  n^)tiuq]dlf  j>  .  (A3) 


Matrix  elements  for  perturbation  operator  Hi  (Appendix  B)  are  given  by 

1 


({Wql^nokwI^/lV’nokoW,!^})  =  \l TT 

V  •  ,=1,2  ^ 


Vk^no  [k(t)j 


+  (qz/qx)2]1,A 

{^q^q})  ^k0,k-q's  A  \jn^  A  ^(^nq/  ,n°,  +li  {^nq,n°})  ^k(),k+q' 

where  a  =  e2 /Tic  is  the  fine  structure  constant. 


x 


,  (A4) 
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For  shortness  of  notations,  Eq.  (A3)  can  be  written  as 

A,(k, «)  -  4  E  /tt'A.Oh »)  =  t) , 


/kk'  =  ^kk'  exP  [£n0(k'(t'))  -  eno(\t(t'))\dt'  j  ,  (A6) 

where  Vkk'  is  defined  in  Eq.  (BIO),  and 

X 

A°(M)  =  T^({^qJ},V'nok(t)|^/|V'nok;(t),«j}) 

xexp  It  [£no(k0(t'))  -  £n0(k(t'))  +  y^(n°i7-  -  nqj)/ia;q]dt,|  .  (A7) 

Also  to  simplify  notations,  we  designated  j(k,f)  =  An(k, t).  Then  separating  terms 
with  diagonal  and  off-diagonal  matrix  elements  /k  k'  in  Eq.  (A5)  so  that 

An(k,  t)  -  ^/kkAn(k,  t)  -  /kk'^n(k',  t)  =  A°n( k,  t) ,  (A8) 

kVk 

the  equation  for  An(k\f)  in  the  third  term  on  the  left-hand  side  of  (A8)  can  be  written  in 
a  similar  fashion 

Ai(k\*)  -  ^/k'k^n(k,t)  -  ^7  /k'k"^n(k",t)  =  A°(k',t) .  (A9) 

’  k'Vk 

Following  to  the  Wigner-Weiskopf-like  approximation,  we  ignore  the  third  term  on  the  right- 
hand  side  of  (A9);  then,  after  integration  on  time,  we  find 


A„(k \t)  =  —  /  dt  fk>k(t)An(k,t)  +  A°(k',t), 


(A10) 


where  we  have  taken  into  account  that  An(k,,t0)  =  <5k'ko,  i-e->  is  eclual  to  zero  for  k  ^  k0, 


<(M)  =  /  A°n(k,t)dt.  (All) 

■h0 

Thus,  by  making  use  of  (A10)  in  Eq.  (A8),  we  obtain  the  latter  equation  in  the  form 

An(k,  t)  -  ^/kkAn(k,  t)  +  ^k'  (*)  [  dtfL'  (*X( k,  t) 


^(k’^  +  ^E-^kk'W^k'.t).  (A12) 
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APPENDIX  B:  MATRIX  ELEMENTS  OF  PERTURBATION  HAMILTONIAN 


Matrix  elements  of  perturbation  Hamiltonian  Hi  for  the  perturbing  quantum  radiation 
field  are  calculated  by  substituting  Hj  given  in  Sec.  II,  with  the  vector  potential  Ar  of 
Eq.  (2),  and  the  wave  function  ipno k(q;  then  we  obtain 


I  Anhc2 


E 


m0cV  u>ceV  [1  +  (qz'/qx)2}1/4 


x 


In  ,  (Sn  , 1’  {^nq,riq})  4  k'  \!  n’ci  i,n  ,+ 1>  {^nq,riq})  4c, k' 

H  q  V  ^  q  q 


(Bl) 


where  the  first  term  on  the  right-hand  side  of  Eq.  (Bl)  is  due  to  absorption  of  a  photon 
from  the  radiation  field  by  the  Bloch  electron  and  the  second  term  represents  emission  of  a 
photon  by  the  Bloch  electron  to  the  radiation  held;  also  we  denoted  (5n  ;  n>  ±1;  {5„q,n'  })  = 

^ n  i,n  ,±  1  riq^q'  4iq,nq  and 


q  q 


q  q 


45  =  /  cos (qxx)  (p  +  Pc)ye±iqzZ^nok'  dr  ,  (B2) 

where  (p  +  pc)y  is  the  y  component  of  [p  +  pc(t)].  In  using  '0nok(t)(r,  t)  = 
n_1/2elK'rMnok(t)(r,  t) ,  the  integrals  /^,  are  equal  to  zero  except  when  K  —  K  ±  qs  =  G, 
where  G  is  the  vector  of  reciprocal  lattice  and  qs  =  {±qx,0,qz}  with  ”  +  ”  for  s  —  1  and 
”  —  ”  for  s  =  2  (Here  we  omit,  for  brevity,  the  prime  on  the  q  vector).  For  values  of  K 
and  qs  in  the  first  Brillouin  zone,  we  can  take  G  =  0.  Then,  it  follows  that  K  =  K  =|=  qs, 
and  since  k (t)  =  K  +  pc(t)/h,  then  we  get  k  =  k  qs.  In  this  case,  the  integrals  1^,  are 
obtained  as 


1 

2 


-4s  ^k',k=Fqs  ’ 

s= 1,2 


(B3) 


with 

4s(±)  =  ^  J (elKrouno k)*(p  +  pc)2/elKro'nnokTqsdr0  ,  (B4) 

where  the  integration  over  dro  is  carried  out  over  the  primitive  cell  volume  Q(>  In  the 
long- wavelength  limit  ( qa  Cl),  the  integrals  zi^  can  be  approximated  as 


l(s±]  =  ^  I  (e*Kr°Mnok)*(p  +  Pc)ye*Kr°nnokdr0  =  J  <ok(p  +  Pc^nok*  • 

In  noting  that  the  well-known  momentum8  expectation  value  for  Bloch  states  is 


(B5) 


(e‘k(',rti„ok)*p  (e‘k<‘,rti„0k)*  =  m0v[k(t)] , 


(B6) 
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where  v[k(f)]  =  (l/ft)VK£n0(K)|k(t),  it  then  follows,  for  k(t)  =  K  +  p C{t)/Ti  in  Eq.  (B6), 
that 

J  Co k(p  +  Pc)  Cok dr  =  m0v(k) .  (B7) 

Thus,  Eqs.  (B3)  -  (B7)  become 

=  ^”Wk)  Y  5k',kTqs  •  (B8) 

s-  1.2 

Then,  using  (Bl)  and  (B8),  the  matrix  elements  of  Hj  are  established  for  use  in  Eq.  (A4). 

Matrix  elements  of  the  perturbing  inhomogeneity  V  are  diagonal  with  respect  to  the 
photon  quantum  numbers 

(M,  C0k(t)ICC0k'(t)’  W)  =  ’  (B9) 

where  we  denoted  for  the  matrix  elements  of  V  between  the  instantaneous  eigenstates 

^kk'  =  (C0k(t),^C0k'(t))  =  [  C0k(qCM)Cok'(t)dr.  (BIO) 
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FIGURES 


FIG.  1:  Dependence  of  the  relative  probability  spectral  density  of  spontaneous  emission  r](ujq) 
[transfer  function  of  Eq.  (23)]  on  the  normalized  photon  frequency  uiq/u jb  calculated  at  resonable 
scattering  tvb  =  10  for  N  =  10. 
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FIG.  2:  Dependence  of  the  relative  probability  spectral  density  of  spontaneous  emission  r}(ujq) 
[transfer  function  of  Eq.  (23)]  on  the  normalized  photon  frequency  lOq/lob  calculated  at  strong 
scattering  tvb  =  2.5  for  N  =  10. 
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FIG.  3:  Dependence  of  the  transfer  function  rjmax  calculated  at  maximum  growth  conditions  of 
Eq.(24)  as  a  function  of  number  of  Bloch  oscillations  N  at  different  values  of  dephasing  parameter 
tub  [5.0  <  tub  <  10.0].  The  dashed  curve  shows  the  dependence  of  r]max  ~  IV2  corresponding  to 
tub  — >  oo,  i.e.,  in  the  absence  of  dephasing  effects. 
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FIG.  4:  Dependence  of  the  transfer  function  rjmax  calculated  at  maximum  growth  conditions  of 
Eq.(24)  as  a  function  of  number  of  Bloch  oscillations  N  at  different  values  of  dephasing  parameter 
tub  [1.0  <  tub  <  4.0].  The  dashed  curve  shows  the  dependence  of  r)max  ~  N2  corresponding  to 
tub  — >  oo,  i.e.,  in  the  absence  of  dephasing  effects. 
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FIG.  5:  Dependence  of  the  prefactor  In(tub)  in  expression  for  the  relative  spontaneous  emission 
probability  at  fundamental  Bloch  frequency  into  the  waveguide  TE\ o  mode  PJe  [Eq. (42)]  calculated 
as  a  function  of  number  of  Bloch  oscillations  N  at  different  values  of  dephasing  parameter  tub- 
The  straight  line  shows  the  dependence  of  In{tub)  ~  N  corresponding  to  tub  — >  oo,  i.e.,  in  the 
absence  of  dephasing  effects. 
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FIG.  6:  Dependence  of  the  relative  spontaneous  emission  probability  PJe  [Eq. (42)]  at  fundamen¬ 
tal  Bloch  frequency  into  the  waveguide  TE\q  mode  calculated  as  a  function  of  number  of  Bloch 
oscillations  N  at  different  values  of  dephasing  parameter  tub-  -PJe  =  1  ( N  >  1)  corresponds  to 
tub  —>  oo,  i.e.,  the  absence  of  dephasing  effects. 


27 


PHYSICAL  REVIEW  B  75,  045330  (2007) 


Microcavity  enhancement  of  spontaneous  emission  for  Bloch  oscillations 

V.  N.  Sokolov*  and  G.  J.  Iafrate 

Department  of  Electrical  and  Computer  Engineering,  North  Carolina  State  University,  Raleigh,  North  Carolina  27695-8617,  USA 

J.  B.  Krieger 

Department  of  Physics,  Brooklyn  College,  CUNY,  Brooklyn,  New  York  11210,  USA 
(Received  15  August  2006;  revised  manuscript  received  2  November  2006;  published  19  January  2007) 

A  theory  for  the  spontaneous  emission  of  a  Bloch  electron  traversing  a  single  energy  miniband  of  a  super¬ 
lattice  while  accelerating  under  the  influence  of  a  constant  external  electric  held  and  radiating  into  a  micro¬ 
cavity  is  presented.  In  the  analysis,  the  quantum  electromagnetic  radiation  held  is  described  by  the  dominant 
microcavity  TE10  rectangular  waveguide  mode  in  the  Coulomb  gauge,  and  the  instantaneous  eigenstates  of  the 
Bloch  Hamiltonian  are  utilized  as  the  basis  states  in  describing  the  Bloch  electron  dynamics  to  all  orders  in  the 
constant  external  electric  held.  The  results  show  that  the  spontaneous  emission  amplitude,  when  analyzed  over 
many  integral  multiple  values  of  the  Bloch  period,  gives  rise  to  selection  rules  for  photon  emission  in  both 
frequency  and  wave  number  with  preferred  transitions  at  the  Wannier-Stark  ladder  levels.  From  these  selection 
rules,  the  total  spontaneous  emission  probability  is  derived  to  hrst-order  perturbation  theory  in  the  quantized 
radiation  held.  It  is  shown  that  the  power  radiated  into  the  dominant  TE10  waveguide  mode  can  be  enhanced 
by  an  order  of  magnitude  over  the  free-space  value  by  tuning  the  Bloch  frequency  to  align  with  the  waveguide 
spectral  density  peak.  A  general  expression  for  the  total  spontaneous  emission  probability  is  obtained  in  terms 
of  arbitrary  superlattice  single  band  parameters,  showing  multiharmonic  behavior  and  cavity  tuning  properties. 

For  GaAs-based  superlattices,  described  in  the  nearest-neighbor  tight-binding  approximation,  the  power  radi¬ 
ated  into  the  waveguide  from  spontaneous  emission  due  to  Bloch  oscillations  in  the  terahertz  frequency  range 
is  estimated  to  be  several  microwatts. 

DOI:  10. 1103/PhysRevB. 75. 045330  PACS  number(s):  73.63.Hs,  72.10.Bg,  73.21.Cd,  73.50.Mx 


I.  INTRODUCTION 

A  theory  of  spontaneous  emission  (SE),  including  the 
probability  for  absorption,  for  a  Bloch  electron  traversing  a 
single  energy  band  in  an  external  constant  electric  held  has 
been  examined  recently  by  the  authors  and  colleagues.1  The 
theoretical  analysis  was  fully  quantum-mechanical,  in  that 
the  radiation  field  was  described  as  the,  free-space  quantized 
electromagnetic  field  in  the  Coulomb  gauge.  Analysis  of  the 
probability  amplitude,  over  integral  multiples  of  the  Bloch 
period,  resulted  in  selection  rules  for  photon  emission  in 
both  photon  frequency  and  wave  vector,  showing  preferred 
transitions  to  the  Wannier-Stark  ladder  levels.  Using  these 
selection  rules,  total  SE  probability  was  derived  to  first-order 
perturbation  theory  in  the  quantized  radiation  field.  Although 
the  output  frequency  of  the  radiation  could  be  operationally 
tuned  to  span  the  gigahertz  to  terahertz  spectral  range  by 
appropriately  fixing  the  constant  electric  held,  the  power  out¬ 
put  for  terahertz  emission  into  free  space  for  a  GaAs-based 
superlattice  (SL)  was  estimated  to  be  about  one-tenth  of  a 
microwatt. 

In  this  paper,  we  analyze  the  SE  of  radiation  emitted  by  a 
Bloch  electron  accelerated  through  the  single  miniband  of  a 
SL  structure,  but  now  with  the  SL  placed  in  a  resonant  mi¬ 
crocavity.  The  intent  of  the  resonant  microcavity  is  well 
known,2,3  to  redistribute  the  free-space  modal  spectral  den¬ 
sity  so  as  to  increase  its  value  at  some  frequencies  and  to 
decrease  it  at  others;  therefore  for  an  active  medium,  such  as 
a  radiating  SL  placed  in  a  cavity  structure,  the  SE  rate  can  be 
enhanced  or  diminished  depending  upon  the  tuning  of  the 
emission  frequency  relative  to  the  cavity  mode  spectral  den¬ 
sity  peak.  We  use  this  tuning  property  offered  by  microcavi¬ 


ties  to  increase  the  power  output  for  enhanced  SE  from 
GaAs-based  SLs.  As  examples  relevant  to  the  principle  idea 
of  this  work,  many  efforts  have  been  focused  on  increasing 
of  the  SE  rate  in  optical  microcavities.4  An  increase  of  the 
overall  emitted  THz  power  of  more  than  one  order  of  mag¬ 
nitude  has  been  reported  by  placing  a  surface-field  emitter 
inside  a  THz  cavity.5 

II.  QUANTUM  DYNAMICS  AND  THEORETICAL 
APPROACH 

For  purposes  of  the  calculation,  we  assume  that  the  SL 
structure  is  placed  in  a  waveguide  with  rectangular  cross 
section  Lx  X  Ly  and  length  Lz,  where  the  coordinate  axes  are 
chosen  to  be  along  the  waveguide  edges.  The  constant  or  dc 
electric  field  E  is  applied  along  the  y  axis,  which  is  also  the 
SL  growth  direction.  The  electromagnetic  field  inside  the 
waveguide,  with  assumed  perfectly  conducting  walls,  is  de¬ 
termined  by  the  guided  modes  corresponding  to  standing 
waves  with  respect  to  the  X  and  Y  axes  [designated  by  an 
integer  pair  (m,n)\,  and  propagating  waves  along  the  Z  axis 
characterized  by  propagation  constant  qz.  Such  modes  form  a 
complete  and  orthogonal  basis  set  for  describing  the  electro¬ 
magnetic  field  within  the  waveguide.  In  the  following,  we 
will  consider  only  transverse  electric  (TE)  modes,  where  the 
electric  field  is  perpendicular  to  the  direction  of  propagation. 
For  practical  cases,6  the  most  important  of  all  confined 
modes  in  this  waveguide  configuration  is  the  TE10  mode 
(m=  1,  n  =  0),  which  is  the  dominant  mode  of  a  waveguide 
with  Lx>Ly.  This  mode  gives  rise  to  the  lowest  attenuation, 
and  the  corresponding  electric  field  E,.,  for  the  chosen  system 
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geometry,  is  polarized  in  the  direction  of  the  dc  held  E. 
Therefore  we  consider  only  one  excited,  TE10,  mode,  and 
ignore  all  the  other  less  effective  TE  and  TM  modes.  For  the 
remaining  TE10  mode,  the  following  conditions  on  E,.  and  H, 
are  in  effect,  namely  that  Erx=Er  z= 0  and  magnetic  held 
component  Hr  =0.  Therefore  it  follows  that  the  vector  po¬ 
tential  A,,  for  the  waveguide  held  has  only  one  nonzero  (y) 
component,  namely,7 


V4  jy  ^ 

- —  sin(  qpc)  {aqe,q-z  +  die ~,q^z) ,  ( 1 ) 

wqey 

where  d 7  and  dq  are  the  photon  boson  creation  and  annihila¬ 
tion  operators,  qx=Tr/Lx,  c  is  the  velocity  of  light  in  vacuum, 
V=LxLyLv  and  e  is  the  dielectric  constant  of  the  medium 
filling  the  waveguide.  The  normalization  constant  in  Eq.  (1) 
is  chosen  in  such  a  way  that  the  Hamiltonian  for  the  quan¬ 
tized  radiation  held  has  the  form  Hr='2qha>qdqdq,  where 
u>q=  wc[l  +  (g'j;/g,J2]1/2  is  the  mode  dispersion  relation,  and 
a>c=qxc/ \s  is  the  angular  cutoff  frequency  determined  by 
the  waveguide  geometry.  The  guided  mode  wavelength  is 
written  as  X  =  X.c/[(<w9/wc.)2- 1]1/2,  where  \C=2LX  is  the  cut¬ 
off  wavelength.8 

The  Bloch  dynamical  properties  are  now  considered  for 
the  situation  in  which  the  electron  is  conhned  to  a  single 
miniband,  “«0,”  of  a  SL  with  band  energy  e„o(K),  while  the 
effects  of  interband  coupling 9  and  electron  intraband  scatter¬ 
ing  are  ignored.  Therefore  the  quantum  dynamics  is  de¬ 
scribed  by  the  time-dependent  Schrodinger  equation 


l^«0W>  =  2  2  AK/Ml^n0k(t),{«qj}> 


k 


X  exp-^ 


if 


e„0(k(?'))  +  2  he 
q  J 


dt'y 

(3) 


where  the  summation  over  k  is  carried  out  over  the  entire 
Brillouin  zone,  and  {«q  ;}  is  specified  over  all  possible  com¬ 
binations  of  photon  occupation  number  nqj  with  photon 
wave  vectors  q  and  certain  polarization  (J=  1,2).  The  instan¬ 
taneous  eigenstates  of  H0  are  given  by9  tpn  k(,)(r,f) 
=  n_1/2e'KrM„ok(f)(r,?),  where  M„ok(f)(r,r)  is  the  periodic  part 
of  the  Bloch  function,  k(f)  =  K+pf.(f)//i,  and  the  values  of 
the  electron  wave  vector  K  are  determined  by  the  periodic 
boundary  conditions  of  the  periodic  crystal  of  volume  ll. 


III.  PROBABILITY  AMPLITUDES,  SELECTION  RULES, 
AND  TOTAL  SPONTANEOUS  EMISSION 
PROBABILITY 

A.  Probability  amplitudes — one-photon  emission 

For  the  case  of  one-photon  SE,  which  assumes  that  ini¬ 
tially  no  photons  are  present  in  the  radiation  field,  the  prob¬ 
ability  amplitude  in  the  wave  function  of  Eq.  (3)  satisfies  the 
initial  condition  A{„  y(k ,  r0)  ={  <5„  ,o}^k,k0  at  time  t=to  when 
the  external  electric  field  is  turned  on.  Here,  K0  and  n°;=0 
are  the  initial  values  of  K  and  nqr  The  probability  amplitude 
for  SE,  A?(k0,r),  at  any  time  t,  is  now  evaluated  in  first-order 
perturbation  theory1  as 


ihjt\Vno(t))  =  H\Vno(t)),  (2) 

where  the  exact  Hamiltonian  H=  [p  -  {el c)  A]2 / 2 m0  +  yr(r) 
+Hr  can  be  reduced  to  a  sum  of  the  following  separate 
Hamiltonians  with  H=H0+Hr+H/.]  Here  the  first  two  terms 
represent  the  Hamiltonian  //0(r)  =  [p+pc(?)]2/2m0+yc.(r),  for 
a  single  electron  in  a  periodic  crystal  potential,  Vc(r),  inter¬ 
acting  with  a  homogeneous  electric  field,  and  the  Hamil¬ 
tonian  Hr  for  the  cavity  mode  electromagnetic  radiation  field. 
The  total  vector  potential  in  the  exact  Hamiltonian  consists 
of  A=AC+Ar,  where  At,=-(c/e)pt,  describes  the  external 
electric  field  with  \>c(t)  =  ef,lE*{t')dt\  and  A,.,  given  in  Eq. 
(1),  describes  the  cavity  mode  quantized  radiation  field;  also, 
m0  is  the  free-electron  mass.  The  Hamiltonian 
H7(r)  =  -(e/m0c)A,.-[p  +  pc(f)],  for  the  first-order  interaction 
between  the  quantum  field  and  the  Bloch  electron,  couples 
both  subsystems  H0  and  and  causes  transitions  between 
the  accelerated  Bloch  electron  states  through  photon  absorp¬ 
tion  and  emission.  Then,  starting  with  the  reduced  Hamil¬ 
tonian  H=Hu+Hr+HI,  use  is  made  of  first-order  time- 
dependent  perturbation  theory  to  calculate  SE  transitions 
probabilities  between  states  of  H0+Hr  while  regarding 
Hj{t)  ~  A,. -[p  +  pc(t)]  as  a  perturbation.10  The  solution  to 
I'P,,  (f))  of  Eq.  (2)  can  be  represented  in  terms  of  eigenstates 
of  basis  states  |  <A„()k(rj ,  {nqJ})  =  |  <A„0k(;)>  I  {«qj})  of  the  unper¬ 
turbed  Hamiltonian  H0+Hr  as 


A3(k0,f)  =  D{qx/q)112  f  dt'vy( k0  -  qs) 

J'o 

[s„0(k0)  -  e„o(k0  -  qs) 

-  ftwjdnj,  (4) 

where  D  =  -ffrrrcal u>ceV,  a=e2! he  is  the  fine-structure  con¬ 
stant,  k0(f)  =  K0+pc(f)//I,  qJ={±<7v,0,^;.}  with  the  “+”  used 
for  ,v=l  and  the  used  for  s= 2,  uv(k(r)) 

=  (1/  h)VK  e„  (K)|k(,p  the  y  component  of  Bloch  velocity  in 
the  band,  and  q=(q2+q2)l/2.  Then  the  spontaneous  emission 
process  results  in  the  total  SE  probability  at  a  time  t,  Pse(t), 
given  by 

pse(t)  =  2  2  K(k0,t)|2.  (5) 

q  5=1,2 

In  evaluating  A?(k0,f),  we  take  into  account  that  the  ex¬ 
ternal  dc  field  E  is  along  the  Y  axis;  then,  it  follows  that 
k0v(t)  =  K0y  +  eE(t-t0)lh.  In  taking  advantage  of  the  periodic 
properties  of  the  terms  in  Eq.  (4),  A?(k0,f)  is  easily  evaluated 
in  terms  of  clocked  integral  multiples  of  the  Bloch  period; 
here,  t=NrB,  where  tb  =  2tt/ wB,  the  time  to  traverse  one  pe¬ 
riod  of  the  Brillouin  zone.  Then  the  time  integral  in  Eq.  (4) 
can  be  replaced  by  an  integral  over  k0v  through  the  substitu- 
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tion  dt=(fil eE)dk0y,  and  the  probability  amplitude,  at  inte¬ 
gral  multiples  of  the  Bloch  period,  can  be  expressed  in  terms 
of  an  integral  over  the  single  Bloch  period,  tb.1’9  Thus  we 
obtain 

\Aq(k0,NTB)\2  =  Sm  2 /l?  lAg(ko» Tb)|2.  (6) 

Sill  \PqlZ) 

where  the  parameter  f3q  is  given  by 

W(  1  fK0y+Gy 

Pq  =  2-77^  +  —  dk0y[sn  Qi0)  -  e  (k0  -  q,)] , 

“b  eEjK0y 

(7) 

and  Gv  =  27r/a,  the  y  component  of  the  SL  reciprocal-lattice 
vector. 


B.  Selection  rules 


From  Eq.  (6),  it  is  seen  that  the  quantity  |A9(k0,A^Tg)|2 
will  reach  its  maximum  growth  value  when  (3q  =  2Tr(m  +  S), 
where  m  is  a  nonzero  integer  and  <5— >0;  for  this  limit,  the 
function  rj((oq)  =  sm2(N/3q/2)/sm2(f3q/2)  — tN2,  i.e.,  it  be¬ 
comes  sharply  peaked  at  the  resonances  (3q=2mn  with  in¬ 
creasing  N.  It  is  clear  that  this  condition  for  maximum 
growth  establishes  the  selection  rule}’9  for  both  the  photon 
emission  frequency,  u>q,  and  the  key  wave-vector  component, 
q7.  Indeed,  from  the  condition  f3q=2,nm,  it  follows  from  Eq. 
(7)  in  the  radiative  long- wavelength  limit  {qa- C  1)  that  one 
generally11  has 


(oq  =  m(oB,  qz=±qzm. 


this  gives  two  conditions,  the  “ Wannier-Stark  ladder”  reso¬ 
nance  frequency  condition,  and  the  sustaining  wave-vector 
cavity  resonance  condition.  These  quantization  conditions 
are  obtained,  naturally,  through  the  use  of  instantaneous 
eigenstates  of  the  Bloch  Hamiltonian  without  requiring  any 
ad  hoc  assumptions  concerning  the  existence  of  Wannier- 
Stark  energy  states.  It  is  worth  emphasizing  that  Eq.  (6) 
should  be  considered  simultaneously  with  Eq.  (8),  i.e.,  both 
selection  rule  conditions  for  the  transition  of  photon  fre¬ 
quency  and  wave  vector  should  be  taken  into  account  explic¬ 
itly  in  evaluating  the  total  SE  probability  spectra  given  by 
|A?(k0,  tb)\2.  In  this  regard,  we  note  that  although  Eq.  (6) 
predicts  the  equal  relative  probabilities  calculated  at  times  t 
=  tb  and  t=NrB  for  different  harmonics  of  Bloch  frequency 
through 


lA,(k0,jVTfi)l2 

K(k0,rB)|2 


=  97  (coq  =  ffldij)  =  N2, 


(9) 


it  is  also  clear  from  Eqs.  (4)  and  (8)  explicitly  that  the  prob¬ 
ability  amplitude  depends  on  a>B  and  q  independently.  There¬ 
fore  the  total  SE  probability  in  Eq.  (5)  must  incorporate,  and 
will  reflect,  the  importance  of  both  selection  rules. 

In  summary,  it  is  noted  here  that  the  probability  amplitude 
of  Eq.  (4)  is  markedly  influenced  by  the  time-dependent  ex¬ 
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ponential  phase  factor.  This  phase  factor  measures  the  elec¬ 
tron’s  field-dependent  energy  change  from  state  e„o[k0(f)]  to 
state  £„  [k0(f)-qs],  while  simultaneously  emitting  a  single 
photon  hojq  to  the  initially  prepared  vacuum  state  of  the  ra¬ 
diation  field.  This  form  of  the  probability  amplitude  does  not 
follow  the  usual  golden  rule  for  stationary  initial-final  state 
dependence  because  the  instantaneous  Bloch  eigenstates  as 
well  as  the  Hamiltonian  describing  the  external  electric  field 
based  on  the  vector  potential  gauge  are  explicitly  time  de¬ 
pendent.  Therefore  as  an  equivalent  alternative  to  the  golden 
rule,  use  is  made  of  the  periodic  nature  of  the  energy  band 
under  consideration,  and  the  maximum  growth  condition  for 
the  probability  amplitude  is  established  at  integral  multiples 
of  the  Bloch  period,  Ntb  with  N  large.  Thus  maximum 
growth  in  probability  amplitude  places  a  quantization  condi¬ 
tion  on  f3q  in  Eq.  (7),  namely,  that  f3q  =  2mn,  which  results  in 
the  selection  rules  of  Eq.  (8)  for  u>q  and  q,,  and  allows  for  the 
calculation  of  the  total  spontaneous  emission  probability. 


C.  Total  spontaneous  emission  probability 

The  total  SE  probability  is  evaluated  at  time  t=NrB,  Pse 
=  Pse(NTB),  by  substituting  |A9(k0,ArB)|2  from  Eq.  (6)  into 
Eq.  (5),  and  then  summing  over  q.  The  sum  over  q  in  Eq.  (5) 
has  been  replaced  by  an  integral  over  q,  taking  into  account 
the  TE10  mode  density  of  states  and  polarization  such  that 
^q(' tt)J  dq(- ■  ■)/[l-(qx/q)2]V2.  The  integral  can 
be  evaluated  by  using  the  property  of  the  integrand  which 
contains  a  sharply  peaked,  symmetric  function  of  q  at  q 
=qm  =  (q2+q2)il2  [see  Eq.  (8)].  Thus  at  every  node  defined 
by  the  resonance  conditions,  the  slowly  varying  function  of  q 
in  the  integrand  can  be  replaced  by  its  value  evaluated  at  q 
-qm,  and  then  removed  from  the  integral  over  q:  after  that, 
the  remaining  integral  can  be  evaluated  to  obtain 


K,(k0  ’Tb)\2 
LxcncHd'-(dM2Y12 


(10) 


Here  lmax  follows  from  qmax = lmax( u>H /  a>c) qx,  and  determines 
the  upper  limit  in  the  sum  over  higher  Bloch  oscillation  har¬ 
monics. 

The  calculation  of  Pse  in  Eq.  (10)  requires  the  use  of 
A?(k0,f)  in  Eq.  (4),  evaluated  when  t  =  t0+rB}2  and  at  the 
selection  rules  of  Eq.  (8),  that  is,  when  hutq=mha)B  and  q 
=qm.  In  addition,  the  dependence  upon  q  in  Eq.  (4)  is  made 
explicit  by  invoking  the  assumption  of  photon  long- 
wavelength  limit,  which  is  valid  for  all  periodic  potentials  of 
interest,  even  SLs,  where  q^irla.  Thus  it  follows  from  Eq. 
(4)  that 


K;(k0,Tfl)l2=  q ;  Adi,  l2,  (11) 

'  qi((oB/2TT) 

where 

h  =7-  [  dfyvyifyjexpi-i/fy)  (12) 

2  ttJ_v 

is  the  Zth  Fourier  component  of  the  vfft/},  the  y  component 
of  electron  band  velocity,  and  d k=k0ya',  D  is  defined  with  Eq. 
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(4).  Since  the  electron  velocity  component  uy(#i)  depends 
on  details  of  the  SL  miniband  structure,  it  is  then  clear  that 
the  band-structure  dependence  appears  explicitly  in  the  SE 
probability  amplitude  through  the  integral  /,  of  Eqs.  (11)  and 
(12).  From  this  general  expression  of  Eq.  (11),  it  is  immedi¬ 
ately  apparent  that  for  electron  dynamics  in  a  purely  har¬ 
monic  miniband,  that  is,  one  in  which  only  one  single  Fou¬ 
rier  component  of  velocity,  uy(#^),  is  nonvanishing  so  that 
1=1,  then  only  the  probability  amplitude  corresponding  to 
the  fundamental  Bloch  frequency  contributes  to  the  total 
emission  process. 


IV.  TOTAL  SPONTANEOUS  EMISSION  RADIATION  AND 
BAND  STRUCTURE  CHARACTERISTICS 

The  analysis  for  total  SE  radiation  characteristics  is  now 
developed  by  considering  a  general  form  of  the  electron  en¬ 
ergy  miniband  dispersion  relation  expressed  as 

oo 

e„  (K)  =  e  (0)  +  2  A,sin2  +  e±(K±),  (13) 

i=i  1 


where  e„o(0)  is  the  miniband  edge,  A;  is  the  width  of  the  /th 
miniband  harmonic  of  the  SL,  and  e±(K±)  is  the  contribu¬ 
tion  from  the  perpendicular  components  of  the  band.  Such  a 
form  of  the  energy-band  dispersion  in  the  SL  growth  direc¬ 
tion  generally  includes  long-range  coupling  over  the  neigh¬ 
boring  QWs  with  a  relative  strength  measured  by  the  specific 
value  of  the  ratio  A/+j /A/<  1 ,  which  is  strongly  dependent 
upon  the  extent  of  wave-function  overlap.  In  particular,  for 
the  well-known  case  of  nearest-neighbor  tight-binding 
(NNTB)  approximation,  only  A  |  with  /=  1  with  purely  har¬ 
monic  energy  dispersion  is  considered  significant,  so  that 
next-nearest-neighbor  and  longer  range  QW  wave-function 
overlaps  are  assumed  to  be  negligibly  small.  The  electron 
group  velocity  in  the  general  miniband  of  Eq.  (13),  for  the 
given  Ky  in  the  y  direction,  is  then  given  by  vy{Ky)  =  (l/h) 
X[<9e„  (Ky)/ dKy\=1,™=lVism(laK^,  where  u/=/aA//2^.,  the 
maximum  velocity  associated  with  the  /th  miniband  of  band 
width,  A?.  Substituting  the  expression  for  vy(Ky)  into  Eq. 
(12),  one  can  find  that  //=-/u;/2.  Then  the  probability  am¬ 
plitude  in  Eq.  (11)  reduces  to 

K,(k0,rB)|2  =  rr2|£)|2^y,  (14) 

a>Bl 


where  account  has  been  taken  for  the  wave  vector  q / 
=qxlu>Bl wc.  Then,  from  Eqs.  (10)  and  (14),  again  using  qt 
=qxlu>B! a>c,  the  total  SE  probability  becomes 


Pse  =  2aNem—^^J 


(v,/c)2 


L  wBl=l  /[I  -  (wJIwb)2] 


2-11/2  ■ 


(15) 


This  is  a  general  expression  for  total  SE  probability  which 
contains  contributions  from  higher  harmonics  of  Bloch  fre¬ 
quency.  Note  that  in  the  NNTB  approximation  of  Eq.  (13), 
obtained  by  letting  A/=A151/,  so  that  v  ,=v ,  Su,  it  follows 
from  Eqs.  (15)  that 


Pse(l  =  1)  =  2  aN 


Lxv\  s 1/2  to2 

LyC1  to2(l  -  w2/w2)1/2’ 


(16) 


where  v l=a\l/2h.  For  a  general  SL  miniband,  higher-order 
harmonic  terms  beyond  the  first  harmonic  term  of  the  Bloch 
frequency  appear  in  the  SE  probability  amplitude  [Eq.  (14)], 
and  thus  appear  in  the  total  SE  probability  [Eq.  (15)]  due  to 
a  nonzero  contribution  from  higher  Fourier  components  of 
Vy(-&k)  [Eq.  (12)]  that  are  in  resonance  with  corresponding 
harmonics  of  the  miniband  energy  spectrum.  In  Eq.  (15),  V/ 
=  laAt/2h  indexes  the  strength  of  the  /th  harmonic  through 
the  parameter  Ah  and  /th  harmonic  is  indexed  through  hoH. 
In  general,  high-order  contributions  are  much  weaker  than 
the  fundamental  one  because  the  coefficients  A,  in  Eq.  (13) 
rapidly  decrease  with  increasing  /. 13,14  Note  that  the  SE  prob¬ 
ability  of  the  higher  harmonics  in  Eq.  (15)  can  be  enhanced 
by  tuning  the  emission  frequency  u>q  =  lu>B  to  align  with  the 
spectral  peak  of  the  waveguide  TE10  mode,  so  that  it  is  close 
to  the  waveguide  cutoff  frequency  coc\  then,  the  resonance  at 
the  fundamental  Bloch  frequency  will  be  suppressed  because 
it  is  less  than  the  waveguide  cutoff  frequency  ( a>B<a>c ).  In 
particular,  for  1=2,  the  SE  probability  for  the  second  har¬ 
monic  generation  is  expressed  as 


Pse(l  =  2)  =  aN 


T  2 

LvV7 


e 1/2 ft)2 


LyC  - (ft)c/2ft)s)2] 


2-11/2’ 


(17) 


where  u2  =  aA2lh.  It  is  noted  therefore  that  without  scattering 
higher-order  harmonic  generation  in  the  SE  spectrum  of  the 
accelerated  electron  from  the  SL  miniband  is  a  signature  for 
band  anharmonicity  in  the  SL  dispersion. 

In  noting  that  the  SE  probability  of  Bloch  radiation 
into  free  space  is  given  by  the  expression  P^s 
=  (2t7/3)q,A/(ui/c)2,1  we  can  compare  both  the  probabilities 
for  SE  at  fundamental  Bloch  frequency  into  free  space  and  a 
microcavity  analyzing  the  ratio 


K  _  3LA-ft)2s1/2 

P%  7rLvft)2[l-(ft)2/ft)2)]1/2’ 


(18) 


where  use  has  been  made  of  Eq.  (16)  for  the  Pse=Pse(l=  1). 
The  enhancement  factor  rj  given  in  Eq.  (18)  is  a  function  of 
the  frequency  ratio,  a>B/wc,  and  increases  monotonically  with 
decreasing  ratio;  of  course,  a>B/ a>c  reflects  the  detuning  of 
the  Bloch  frequency  relative  to  the  TE10  mode  cutoff  fre¬ 
quency.  Taking  for  an  estimation  e=12.2  and  Lx/Ly=2,  one 
can  find  that  for  a>B/wc<  1.15  the  enhancement  factor  77  in¬ 
creases  precipitously  with  decreasing  ratio;  at  cjb/wc—  1.15, 
77—  10,  indicating  that  when  wB  is  tuned  to  within  1 . 1 5ftJr, 
the  SE  enhancement  will  be  increased  at  least  by  an  order  of 
magnitude  over  unity. 


V.  SPONTANEOUS  EMISSION  POWER  ESTIMATE 

For  numerical  estimations,  we  assume  a  GaAs-based  SL 
structure  with  the  SL  lattice  parameter  a  =100  A,  vertical 
dimension  9  /urn,  and  lateral  cross  section  18X  1000  /2m2 
embedded  into  a  rectangular  waveguide  with  horizontal  and 
vertical  dimensions  Lx/Ly=2.  The  electron  density  in  the  ac¬ 
tive  region  is  taken  to  be  5  X  1016  cm-3.  Taking  for  the  SL 
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lowest  miniband  energy  width  A  =  20  meV,  the  maximum 
group  velocity  in  the  miniband  is  estimated  as  vmax=l.6 
X  107  cm/s.  These  parameter  magnitudes  are  close  to  those 
of  GaAs-based  SL  structures  used  to  study  high-frequency 
microwave  generation. 15-1 7  Spontaneous  emission  of  a  pho¬ 
ton  with  the  energy  10  meV  corresponds  to  the  Bloch  fre¬ 
quency  vb  =  mbI2tt=2.5  THz.  The  electric  field  required  to 
achieve  such  Bloch  frequency  is  E=  fia>B/ea  =  10  kV/cm, 
and  results  in  the  application  of  9  V  across  the  vertical  di¬ 
mension  of  the  SL  structure.18  The  spontaneous  emission 
probability  of  radiation  into  free  space  can  be  estimated  tak¬ 
ing,  for  example,  N=  100  as  Pfs= 4.3  X  10-7;  and  the  genera¬ 
tion  energy  per  electron  ha>BPSjs=4.3  X  KL6  meV.  Since 
there  is  a  total  of  «  =  8  X  109  electrons  in  the  active  region  of 
the  SL,  the  generated  energy  achievable  is  estimated  to  be 
Pfs  =  nhu>BPjs=34A  eV,  which  corresponds  to  a  power  out¬ 
put  generated  into  free  space  WfS=(vB/N)Pfs—  0.14  /j,W.  Al¬ 
though  the  power  generated  into  free  space  is  discernibly  low 
for  the  SE  of  Bloch  oscillation  radiation,  it  follows  from  Eq. 
(18)  that  SE  probabilities  and  rates  are  substantially  modified 
for  the  case  of  radiation  into  the  waveguide  mode.  We  find 
from  Eq.  (18)  that  77=  20,  if  we  take  for  the  detuning  param¬ 
eter  wB/ coc=  1.05.  Then,  using  the  obtained  value  for  the  en¬ 
hancement  factor,  we  estimate  the  power  output  generated 
into  the  TE10  waveguide  mode  as  Wwg  —  3  //W .  For  purposes 
of  this  estimate,  the  cavity  was  assumed  to  be  lossless;  fur¬ 
ther  studies  will  consider  cavities  with  finite  quality  factor. 

It  is  noted  that  a  Bloch  oscillation  SL  does  not  require 
controlled  inversion  population  between  Wannier-Stark  lad¬ 
der  levels  to  get  the  desired  SE  photon  frequency;  the  desired 
frequency  is  controlled  by  the  applied  field.  Whereas  in  other 
SL  light  generating  devices,  such  as  quantum  cascade  lasers, 
a  large  inversion  population  is  required  to  provide  stimulated 
emission  with  resulting  high  threshold  current  densities  and 
high  heat  dissipation.  In  this  regard,  the  Bloch  oscillator  in 
SE  offers  a  novel  option  for  operating  at  THz  frequencies, 
provided  the  power  output  can  be  enhanced  in  the  coherent 
Bloch  regime. 

VI.  SUMMARY 

A  quantum  theory  has  been  presented  to  describe  sponta¬ 
neous  emission  of  radiation  for  a  Bloch  electron  accelerating 
across  a  single  energy  band  of  a  SL  in  a  constant  electric 
field  while  the  radiating  SL  interacts  with  a  resonant  micro¬ 


cavity.  In  the  theoretical  analysis,  use  was  made  of  the  in¬ 
stantaneous  eigenstates  of  the  Bloch  Hamiltonian  in  the  ex¬ 
ternal  electric  field  described  in  terms  of  the  vector  potential 
gauge,  and  the  spontaneous  emission  rates  were  calculated 
using  first-order  time-dependent  perturbation  theory  in  the 
TE10  waveguide  mode  quantized  radiation  field.  The  analysis 
of  spontaneous  emission  amplitude  resulted  in  transition  se¬ 
lection  rules  for  both  photon  emission  frequency  and  wave 
number  at  the  Wannier-Stark  ladder  levels;  this  result  is 
found  naturally  through  the  use  of  the  instantaneous  eigen¬ 
states  of  the  Bloch  Hamiltonian  with  no  ad  hoc  assumptions 
employed  concerning  the  existence  of  Wannier-Stark  quan¬ 
tized  energy  levels  within  the  band. 

Using  the  selection  rules,  the  total  spontaneous  emission 
probability  was  calculated  for  a  general  SL  miniband  disper¬ 
sion,  and  was  shown  to  depend  upon  the  fundamental  and 
higher-order  Bloch  frequency  harmonics  as  well  as  the  band- 
structure  parameters;  in  the  limit  of  nearest-neighbor  tight 
binding,  the  results  for  total  spontaneous  emission  reduced  to 
the  expected  sole  dependence  upon  the  fundamental  Bloch 
harmonic. 

Finally,  it  was  shown,  within  the  nearest-neighbor  tight- 
binding  model,  that  the  total  spontaneous  emission  probabil¬ 
ity  in  the  microcavity  is  substantially  enhanced  over  the 
comparable  free-space  value  when  the  Bloch  frequency  is 
tuned  by  means  of  the  external  electric  field  so  as  to  align 
with  the  spectral  peak  of  the  waveguide  mode  density  of 
states.  In  this  regard,  a  theoretical  estimate  showed  that 
about  a  one-order-of-magnitude  enhancement  factor  for 
GaAs-based  SLs  resulted  in  a  power  output  of  about  3  /jl W. 
Here,  it  is  noted  that  the  prime  source  of  the  total  spontane¬ 
ous  emission  enhancement  is  due  to  the  alignment  of  the 
Bloch  frequency  with  the  peak  of  the  spectral  density  of 
states;  in  the  spirit  of  the  Purcell  effect,2  this  will  be  a  gen¬ 
eral  feature  of  any  cavity  considered  for  this  analysis. 

Future  directions  of  this  work  include  studying  the  limit¬ 
ing  effects  of  dephasing  inhomogeneities  such  as  SL  inter¬ 
face  scattering  and  other  competitive  processes20-22  so  as  to 
more  realistically  evaluate  the  optimal  magnitudes  of  power 
output  from  the  spontaneous  emission  of  SL  Bloch  oscilla¬ 
tion  radiation. 
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